function [ systoleVel ] = verifySystoleStart ( systoleVel )
% This function identifies the true start of systole

systoleXVals = systoleVel(:, 1);
systoleYVals = systoleVel(:, 2);

% Determines the index and x value of the maximum velocity peak
[ globalMaxIdx, globalMaxXVal ] = findMaxVel( systoleXVals, systoleYVals );

% Extract the curve from the beginning of systole to peak velocity
[ curveVals ] = extractMaxMinCurve ( systoleVel, 1, globalMaxXVal );

% Identify slope of curve
[ curveVals ] = slopeDiff ( curveVals );

% Remove regions at the start of the curve which do not correlate to the
% start of systole based on slope
[ curveVals ] = getPositiveCurve ( curveVals );

systoleVel = vertcat(curveVals, systoleVel(globalMaxIdx+1:length(systoleVel), :));


end

